Large-order shifted 1/N expansions through the asymptotic iteration method 



T. Barakat * 



Physics Department, King saud University, 
Riyadh 11451, Saudi Arabia 



Abstract 

The perturbation technique within the framework of the asymptotic iter- 
ation method is used to obtain large-order shifted 1/N expansions, where iV 
is the number of spatial dimensions. This method is contrary to the usual 
Rayleigh-Schrodinger perturbation theory, no matrix elements need to be cal- 
culated. The method is applied to the Schrodinger equation and the non- 
polynomial potential V(r) = r 2 + r^^pi\ in three dimensions is discussed as 
an illustrative example. 
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1 Introduction 



The shifted 1/N expansion technique (SLNT) proposed by Sukhatme and Imbo 
[1] is an extremely powerful method of solving Schrodinger equation, and has been 
used extensively to determine the eigenenergies for some important potentials [2-6]. 
The SLNT is an expansion in powers of A -1 / 2 , A = N + 21 — a, where N, £ and a 
stand for the number of dimensions, the angular momentum quantum number, and 
a is a properly chosen shift parameter, respectively The shift parameter a is usually 
chosen so as to improve the convergence of the energy perturbation series, and to 
obtain the correct answer for the harmonic oscillator and the hydrogen atom to all 
orders. 

After expanding the potential-energy function and the centrifugal term in Tay- 
lor series about an appropriate point r one is left with the Hamiltonian operator 
for a harmonic oscillator plus a polynomial perturbation. Then one applies pertur- 
bation theory and obtains the perturbation corrections for the eigenf unctions and 
eigenenergies. 

The Rayleigh- Schrodinger and the logarithmic perturbation schemes (referred 
to as RSPT and LPT respectively) have been used for the calculation of these 
corrections. 

As well known the applications of the above two schemes were restricted by 
serious difficulties. They require considerable computational time and effort, and 
they involve, in general, quite elaborate algebraic manipulations, so it was extremely 
laborious to advance beyond the first four perturbation expansion terms in the 
eigenenergy series [2]. 
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On the other hand, considerable progress was made in the last few years to 
obtain large-order shifted 1/N expansions [7-10]. For example, Maluendes et al. 
[7] reported an approach in which the coefficients of the SLNT of arbitrarily high 
orders could be generated by means of the hypervirial (HV) and Hellmann-Feynman 
theorems (HF), and thereby providing an excellent check for the convergence of the 
method. 

However, the previous authors in their work did not give explicit expressions of 
their algorithm, each order getting progressively much more complicated than the 
previous one, and the derivations were tediously long. Thus, the need arises here to 
have a relatively simple, fast and effective method that provides large-order shifted 
1/N expansions. 

The so-called perturbation technique within the framework of the asymptotic 
iteration method (AIM) [11] have emerged in recent years to be a very useful and 
powerful technique of attack. This method reproduced excellent results for many im- 
portant potentials in relativistic and non-relativistic quantum mechanics. Through 
AIM one can actually obtain all the perturbation corrections to both energy level 
shifts and wavefunctions for all states [12, 13]. These quantities can be calculated 
to any given accuracy, since the generation of successive corrections in the present 
perturbative framework, only requires the solution of simple algebraic equation. 

The method is also applicable in the same form to both the ground state and 
excited bound states without involving tedious calculations which appeared in the 
available perturbation theories. 

Encouraged by its satisfactory performance, we feel tempted to extend AIM, and 
to see this time how the AIM can be used to obtain large-order shifted 1/N expan- 
sions for the three-dimensional Schrodinger equation with any arbitrary spherically 
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symmetric potential V(r) directly without either, using the base eigenfunctions of 
the unperturbed problem, or needing to calculate matrix elements. 

As an illustration, the present technique is applied to the non-polynomial po- 
tential V(r) = r 2 + (Y^2j- This potential appears in several areas of physics. In 
field theory it provides a simple zero- dimensional model possessing a non-polynomial 
Lagrangian [14]. In laser physics it arises out of the Fokker-Planck equation for a 
single-mode laser [15]. 

For this potential only a class of exact analytical solutions for certain parameter 
dependence b = 6(c) were obtained [16, 17]. Hence, it has been a subject of several 
investigations and many authors have studied the one-, two-and three-dimensional 
cases [18, 19]. Roy, Roychoudhury and Roy have shown the supersymetric character 
of this potential and given new solutions using the standard 1/N expansion method 
[20]. 

With this in mind, this paper is organized as follows. In Sec. 2, the formulation of 
SLNT through AIM is outlined to find the eigenenergies for any arbitrary spherically 
symmetric potential. The analytical expressions for AIM are cast in such a way that 
allows the reader to use them without proceeding into their derivation. In Sec. 3, 
we explained how to obtain numerically the eigenenergies for the non-polynomial 
potential V(r) = r 2 + nj^y , and therein we will compare the accuracy of our results 
with those obtained by Roy et al. [20]. Finally, the paper ends with a brief summary 
and concluding remarks on the method and our findings. 

2 Formalism of the asymptotic iteration method for SLNT 

The radial part of the time-independent Schrodinger equation for central-field 
model in terms of the expansion parameter A in N- dimensional space with (h = 



3 



d 2 A 2 2A 4B U 

33 + 73( 1 + -T + Tt) \Xnt{r) = E neXn i(r), (1) 



2m = 1) is 



<ir 2 4r 2 A A 2 
with 

A = N + 2£-a, A=l-N + a, B = (N — a)(N — a — 2)/4. 

SLNT begins with shifting the origin of the coordinate through the definition 

x = K l /\--l\ (2) 
f o 

where r is chosen to minimize the effective potential V e ff(r) = A 2 (^a + ^a^), so 
that 

A 2 = 2r V(r ). (3) 
Expansions about r = r , yield 



<i 2 

+ ^(a^V +2 + + &<? i+ V) } x^(x) = e n ,Xn£(^) (4) 



cix 2 

i=0 



where 



a - ( iy^ + 3 i r Q +1 dl+2y ^ 8.-(-iyl±±A e.-(-lY(i + l)B (5) 
«.-(!) 4 + 2 (^ + 2)!V'(roK +2 ' A_l } 2 A * " ( + 1)B > {5} 

and 

e„ £ = rlg\E nl - V eff (r )), and g = 1/A 1 / 2 . (6) 
The systematic procedure of the AIM begins now by expanding the eigenenergy 



4 



term in powers of g, and rewriting equation (4) in the following form 

d 2 



dx 2 



+ a x + fa + g(atix + fax) + g (a 2 x + (3 2 x + £ 



+g 3 {a 3 x 5 + (3 3 x 3 + frx) + g\a 4 x 6 + p A x A + £ 2 x 2 ) + • • • ] Xn e(i 



' (0) , (1) , 2 (2) . 3 (3) , 4 (4) . 



Xn£{X) 



If we further insert the ansatz 



Xm{x) = e lx2/2 fnz{x) 



(7) 



(8) 



into equation (7), carrying out the mathematics, in this case, the function f n i{x) 
will satisfy a new second-order homogeneous linear differential equation of the form 



f'ni( x ) = x o( x ,g)f' n e( x ) + s (x,g)f n e(x), 



(9) 



where X (x,g) = 2^x, and 



s (x, g) = (a - 7 2 )x 2 + f3 + 7 - e nE + g{aix 3 + fax) + g 2 (a 2 x 4 + f3 2 x 2 + £ ) 
+ g 3 (a 3 x 5 + (3 3 x 3 + £ 1 x)+g 4 (a i x 6 + f3 4 x 4 + £ 2 x 2 ) + --- , (10) 



c - M J_ n J& J_ r?J^ _L „3j3) 4 (4) 
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Here, it should be pointed out that, the choice of g in equation (9) will be 
motivated, that is when we switch off g, equation (9) will be reduced to an exactly 
solvable eigenvalue problem within the framework of AIM [21-23]. 

To apply the perturbation expansion technique within the framework of AIM we 
rely on the symmetric structure of the right hand side of equation (9). Thus, we 



differentiate equation (9) (k + 2) times with respect to x, k — 1, 2, .... Then we take 
the ratio of the (k + 2) th and (k + l) th derivatives, and for sufficiently large k we 
introduce, respectively, the "asymptotic" aspect and the termination condition of 
the method, which, in turn will lead to 

,(.,^£^=^44 (12) 

5 k (x, g) = s k (x, g)\ k+1 (x, g) - s k+1 (x, g)\ k (x, g) = 0. (13) 

We now proceed to obtain the eigenenergies of equation (9) systematically in terms 
of the expansion parameter g. If we expand 8 k (x,g) around g = 0, we get the 
following series 

x ( x x , n x , 9 d5 k (x,g) g 2 d 2 5 k (x,g) g 3 d 3 5 k (x,g) 

S k (x,g) = S k (x, 0) + - — \ g=0 + - — | fl=0 + ^7 dg3 \ g =o + --(14) 

According to the procedure of AIM [21-23], 5 k (x,g) must be zero; if this to be true 
for every g value, then every term of the series must be zero. That is to say 

^ (li9) = J*4fed u = 0i J = 0>1>2 .„. (15) 

A quantitative estimate for e n e expansion terms can be obtained by comparing 
the terms with the same order of g in equations (9) and (14). Therefore, it is clear 
that the roots of 5 k °\x,0) = give us the zero'th contribution energy terms e^j. 
Likewise, the roots of 5 ( k 1 \x, g)\ g=0 = give us the first correction terms e^, and 
so on. Therefore, the general solution for the eigenenergies E n e in conjunction with 
equations (6) and (11) is 



E n o — 
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3 Numerical results for the eigenenergies of the potential V[r) = r 2 + 

Within the framework of the asymptotic iteration method mentioned in the above 
section, the eigenenergies E ne of the non-polynomial potential V(r) = r 2 + 
are calculated by means of equation (16). 

To obtain the zero'th contribution energy terms e^J, one should simply switch 
off g in equation (9), that will lead to an exactly solvable eigenvalue problem within 
the framework of AIM, 

C(x) = Tryxfjx) + ((a„ - 7 V + A) + 7 - fm(x). (17) 

For each iteration, the expression 5^\x,0) = in equation (13) depends on two 
variables namely 7 and x. Since the problem is exactly solvable, the calculated 
eigenenergies e^j by means of this condition are independent of the choice of x once 
we set 7 = y/oo; then the roots of 5^ (x, 0) = are 

ejg ={3 + (2n+l)(a ) 1/2 , n = 0,1,2, .... (18) 

As we noted before, the leading contribution term of the total energy is of order 
A 2 given in equation (16). The next contribution is of order A and is given by 
Po + (2n + l)(a o y/ 2 . It is customary to choose the shift parameter a so as to make 
this contribution vanish. This choice is physically motivated by requiring agreement 
between the 1/A expansions and the exact analytic results for the harmonic-oscillator 
and Coulomb potentials to all orders [2]. However, there are some difficult cases in 
which this simple choice is insufficient and it is really necessary to select an order- 
dependent value of a according to minimal sensitivity or other appropriate criterion 
[7]. Nevertheless, at this point, it is enough for us to compare the AIM results with 
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the SLNT results. Therefore, we choose 

a = 2-2(2n + l)(a ) 1/2 , (19) 

A = N + 2£ — a, and A 2 = 2rQV'(r ). Collecting these terms and carrying out the 
mathematics, one can get 

3 + 2 £ - 2 + 2(2n + l)(a ) 1/2 = (2r V'(r )) 1/2 , (20) 

which is an implicit equation for r . Once r is determined, the leading term 
A 2 V e //(r ) can be calculated numerically. On the other hand, to obtain the higher 
order perturbative expansion terms, first one should go to equation (9) switch on 
g and then replace e n £ with e^j + ge^, and terminate the iterations by imposing 
the condition S^\x,g)= as an approximation to equation (9). The first root of 
the resulting equation gives e^j . Similarly, and very easily one can obtain the other 
perturbative expansion terms. 

Throughout the present calculations, it is observed that the perturbation correc- 
tions of odd orders e^ +1 ^ are vanish for all %. 

In table 1 an explicit list of calculations up to six'th order with different values 
of £, b, and c are given, so that the reader may, if so inclined, reproduce our results. 
In table 2 the results of AIM, together with the standard shifted expansion method 
E n £{l/N) and the exact super-symmetric results are displayed for comparison pur- 
poses. 

In tables we have only considered the eigenenergies for the ground state n — 0. 
This was in order to make a clear comparison between the results of this method 
and the results of [20]. Examination of tables shows that the accuracy of the AIM 
is better than the accuracy of shifted 1/N expansion method, and the predicted 
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eigenenergies E^AIM) are all in excellent agreement with the results of the super- 
symmetric method [20]. 

We have also shown that, it is very easy task to implement the perturbation 
technique within the framework of the AIM without having to be worry about the 
ranges of the couplings in the potential. 

This method is also applicable in the same form to both the ground state, and 
excited bound states without involving tedious calculations which appeared in the 
available perturbation theories. 

As a concluding remark, the present method enable one if it is necessary, to keep 
a as a free parameter up to the end of the calculations. With this choice the results 
can be drastically improved by raising up the perturbative order in the expansion 
to any order, and then one can determine the value of a according to the minimal 
sensitivity method or any other appropriate criterion. 
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Table 1: The calculated values of the coefficients in the energy expansion e nl for the 
non-polynomial potential by means of this work with different values of £, b, and c. 



i 


c 


b 


e (o) 


e (1) 


,(2) 


e (3) 


e (4) 


,(5) 


e (6) 





0.1 


-0.46 








-0.021504914 





0.01472806 





-0.010867 


1 


0.1 


-0.5 








-0.030894587 





0.02042523 





-0.0034827 


2 


0.1 


-0.54 








-0.034627947 





0.01793987 





0.018511 





0.01 


-0.0406 








-0.000018842 





1.707706x10-6 





-1.7717x10-6 


1 


0.01 


-0.041 








-0.000049912 





7.312060x10-6 





-7.6296x10-6 


-1 


0.1 


-0.42 








-0.005539564 





0.002341067 





-0.0018839 
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Table 2: Comparison between selected eigenenergies calculated from the standard 
shifted (1/N) expansion method E n e(l/N) [20], the exact super-symmetric values 
E n e(SUSY) [20], and the eigenenergies E n i{AIM) computed by means of equation 
(16) up to four'th and six'th orders. 



i 


c 


b 


E oe (l/N) 


Eoe(SUSY) 


Elf {AIM) 


E^f {AIM) 





0.1 


-0.46 


2.400520 


2.4 


2.40051591814138 


2.3999024133815471 


1 


0.1 


-0.5 


4.000116 


4.0 


4.00012195111641 


4.0000777082899992 


2 


0.1 


-0.54 


5.599965 


5.6 


5.59998537049538 


5.6000735787721094 





0.01 


-0.0406 


2.939999 


2.94 


2.94000001431155 


2.9399998857627638 


1 


0.01 


-0.041 


4.899974 


4.9 


4.90000002136554 


4.8999999017683162 


-1 


0.1 


-0.041 


0.801177 


0.8 


0.80117658318143 


0.7984848338375296 
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